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Abstract — In this paper, we consider the problem of 
optimum multi-domain real-time beamforming and high- 
precision beam pattern positioning in application to very 
large wideband array antennas, particularly to the Square 
Kilometre Array (SKA) aperture array antenna. We present 
a new structure for wideband space-frequency beamforming 
and beamsteering that maximizes detectability of cosmic 
signals over the array operational frequency range. 

1. Introduction 

The SKA is a multi-purpose radio telescope that will operate 
at metre and centimetere wave-lengths with extremely high 
sensitivity and wide field of view (FoV) and will have 
unprecedented surveying power. The range of key science to 
be tackled by the SKA covers the epoch of re -ionization, 
galaxy evolution, dark energy, cosmic magnetism, strong 
field tests of gravity, gravitational wave detection, 
transients, proto-planetary discs, and the search for extra- 
terrestrial life [1]. 

More specifically, the SKA is to be a radio telescope with 

- the sensitivity to detect and image hydrogen in the early 
universe through its enormous collecting area of about 1 
million square metres. This will make it about 50 times 
more sensitive than the Extended Very Large Array 
(EVLA), and able to reach an rms noise level of 10 nano- 
Jy in an 8 hour integration for a continuum observation, 

- 50% of the collecting area concentrated in the central 5 
km diameter for optimal detection of hydrogen, pulsars, 
and magnetic fields, 

- a fast surveying capability over the whole sky through 
its sensitivity and very large angle FoV. The SKA will be 
1 0000 times faster than the EVLA in surveying the sky, 

- the capability for detailed imaging of compact objects 
like active galactic nuclei through its large physical extent 
of at least 3000 km, 

- a frequency range from 70 MHz to 10 GHz, 

- data transport to the central data processor via very 
wideband (terabit/sec) fibre links. 
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Figure 1 - Artist's impression of the SKA central core. 

The SKA is an interferometer array capable of imaging the 
radio sky at frequencies from 70 MHz to 10 GHz, and 
providing an all-sky monitoring capability at frequencies 
below 1.4 GHz. It covers the frequency band with three 
different kinds of receiving systems [2]-[4]: 

1. Sparse aperture arrays (AA-lo) for 70-450 MHz, 

2. Dense aperture array tiles (AA-hi) in the core of the array 
for all-sky monitoring in the frequency range 0.4-1.4 GHz, 

3. A 15m parabolic dish array, with wideband single-pixel 
feeds for 1.2-10 GHz. 

These three components all make use of the same data 
communications, processing, and software. An artist's 
impression of the array central core is shown in Figure 1. 

The overall structure of the SKA signal processing system is 
shown in Figure 2 and includes the collector systems on the 
left, communications and control network in the centre and 
correlation and processing on the right [4]. 

A key cost driver for the AAs is the highest frequency 
supported, due to each element having an effective area 
which is a function of A 2 , hence, the number of elements 
required for a given sensitivity increases quadratically with 
frequency. 
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Figure 2 - SKA implementation using AAs and single pixel feeds on dishes 



Each AA-hi station is a dense ~50m diameter circular plane 
array consisting of N A ~ 75,000 dual polarisation elements. 
An outline design of the AA system is shown in Figure 3. 
The design consists of four main blocks: 

1. The front-end collectors. Each element of the AA-hi and 
AA-lo is positioned as part of the array design and tightly 
designed with its associated LNA for the lowest noise front- 
end design. This is amplified and passed to the "Tile 
processor" for initial beamforming. 

2. Tile processors perform the first stage of beamforming, 
where ~64 dual polarisation elements for the AA-hi are 
arranged using the most effective mix of RF and digital 
techniques, to form a number of tile beams. The bandwidth 
between the Tile Processors and the Station Processors is a 
key determinant of the performance of the AAs. 

3. Station processors. These bring together the output of all 
the AA tiles. They form the beams for transmission to the 
correlator. The calibration algorithms to form high precision 
station beams will be handled primarily by the station 
processors. 

4. The control processors keep the operation of the station 
coupled to the rest of the SKA. They also monitor the health 
of the arrays, detect non-functioning components and adjust 
the calibration parameters appropriately. 

Again, the signal processing is implemented in a two stage 
structure, the Tile Processors of the first stage perform the 
initial digitisation and beamforming, the number of primary 
beams - N B , their frequencies plus bandwidth and bits per 



sample can be configured at the observation time. The 
constraint is the total data rate of the internal digital links. 

Identical beams from all the tiles in each array are then 
combined in the Station Processors to produce the required 
station beams. To achieve the total FoV needed there are 
many hundreds of individual station beams. Every first 
stage processor has a link to all the Station Processors. The 
Station Processors each link to wide area communication 
fibres directly to the central correlator. 

The signal processing, which will be using integer (fixed- 
point) arithmetic for efficiency in power and silicon, 
requires a capability estimated for an AA-hi array of ~10 
Peta-MACs. It has become clear that this performance level 
is achievable using either a dedicated ASIC solution, or 
massively multi-core parallel integer processors. This latter 
processing solution is very attractive for flexibility and the 
implementation of novel algorithms. 

The main objective of the paper is the design and 
performance optimization of a cost-effective structure of all- 
digital real-time space-frequency beamforming and 
beamsteering for the SKA station dense AA. 

As mentioned above, the main challenges of the design are 
associated with limited computational power of a given 
silicon technology. 

The upper bound of the computational power assessed 
through the product of the data bandwidth (data rate) and 
dynamic range (number of bits per sample) is a key 
determinant of all-digital beamforming performance. 
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Figure 3 - Outline AA station 



According to the SKA research and design studies [2]-[4], 
forthcoming DSP technologies will be capable of handling 
the input data ranged from 4-bit at 2.4 GS/s to 6-bit at 3.4 
GS/s over the 700 MHz frequency extent. 

Unfortunately, these data attributes do not allow the SKA 
antenna processor designers to "digitize" directly the analog 
WB/UWB beamforming structures outlined in [5] due to 
enormously high computational power that would require 
for the all-digital implementation. 

The leading idea of realizable all-digital WB beamforming 
underlain in the SKA station AA design [2] -[4] is to reduce 
the internal and output data rate by making use of parallel 
FFT-based frequency splitting applied directly to the AA 
sensor signals. 

2. Frequency-Space Beamforming 

A generalized schematic of the AA-hi beamformer that has 
been developed by the current stage of the SKA design [2]- 
[4] is shown in Figure 4 where N A is the number of antenna 
array sensors, Nb is the number of spatial beams, (0„,8„) is 
the direction the «-th beam (n=l, ... , N B ) is focused on and 
N F is the number of frequency subbands to be specified by 
requirements of radio-astronomic observations. 

The Frequency-Space Beamformer (FS-BF) depicted in the 
figure is a sequence of the following array signal processors: 
the polyphase FFT-based frequency-domain demultiplexer 



purposed to form N F separate frequency slots for each of N A 
array antenna elements, the space-domain beamformer 
aimed to create Nb controlled spatial beams for each of N F 
frequency subbands, and the calibration processor to correct 
amplitude-phase errors in each of N B X N F resulting beams. 

For obvious reasons, the number of linearly-independent 
partial beams, Nb should not to be larger than the number of 
array antenna elements Na, i.e. N B < Na- At the current stage 
of the SKA design study, 

N A ~ 70, 000, Nb~1200 and N F = 1024 ( 1 ) 

Let us assess the total number of fixed-point real-valued 
multiplications the FS-beamformer requires to process N F 
input successive snapshots in forming N B spatial beams for 
each of N F frequency subbands. This value will serve as a 
benchmark in estimation of a beamforming structure's 
implementation cost. 

The number of multiplications required to perform the 
frequency-domain stage of the FS-beamforming is 

M FS1 = 0{2N A N F log 2 N F } (2) 

At the second stage of processing, when the FS-beamformer 
assembles and focuses the corresponding spatial beams, the 
number of real- valued multiplications it consumes is about 



3. Space-Frequency Beamforming 



sc : 



sc : 



FIR Filters 




FIR Filters 


FFT 




FFT 



<4>A>- 




Linear 
Combiners 



(+A>- 



Linear 
Combiners 



Frequency 
splitting 



Space 
beamforming 

and 
beamsteering 



Cal & Corr 



beam 1 beam N B 

subband 1 



7n b 



Cal & Corr 



beam 1 beam N B 

subband N p 



Obeying the theory of optimum space-time signal 
processing [6], we conceive the SKA station beamformer as 
a wideband spatial filter to be matched to any kind of 
expected radio-astronomy signals through maximizing the 
sample signal-to-noise power ratio (SNR) at the filter 
outputs at each time discrete. 

Thus, in order to make the space-time filter really optimum 
for any radio -astronomic signal, both full-band and 
narrowband, we have to arrange spatial beamforming and 
beamsteering just before frequency demultiplexing. 

A newly designed structure of the all-digital wideband AA- 
hi beamformer that could be named as the Space-Frequency 
Beamformer (SF-BF) is presented in Figure 5. 



Figure 4 - Frequency-Space beamformer 

M FS2 = 0{4N B N A N F ) 



(3) 



Thus, the total number of real-valued multiplications the 
FS-beamformer requires is 



M FS =M FSI +M FS2 = 0{2N A Np{\og2NF + 2N B )} 



(4) 



As follows from the theory of optimum space-time signal 
processing [6], the FS-beamforming scheme that has been 
derived using the "narrowband" approach [2]-[4] is a space- 
time filter well-matched only to the signals whose 
autocorrelation function widths are much longer than the 
maximum time delay of the signals envelopes across the 
array aperture. 



substantially non-optimal (lossy) for a wide variety of 
WB/UWB cosmic signals including short-living transients, 
pulsars fine-grained signatures, etc. 
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Figure 5 - Space-Frequency beamformer 



It needs also to add the following remarks concerning HW 
implementations of the FS-beamformer. 

Short bit-grid size quantised sample sequences are known to 
have the quantisation noise component attributed with the 
irregular multi-peak power spectrum. In order to mitigate 
the quantisation aliasing, we have to use an additional bank 
of Nj x Nf passband FIR-filters at the FFT's outputs. 

The phase-shifting realization appears to be numerically 
expensive as it requires two operations of full dimensional 
complex multiplication. 



From the figure is seen that the SF-beamformer contains N B 
parallel processing arms to form N B individually steered 
wideband spatial beams, with each beam decomposed onto 
N F frequency subbands. 

Each beamforming arm of the scheme is a two -stage real- 
time digital signal processor containing the wideband spatial 
processor to form the full-band beam pattern directed 
towards the prescribed point, (ip„,d„), and the polyphase 
FFT-based demultiplexer that performs final frequency 
multiple down-conversion. 



At last, the FS-BF structure suffers from WB frequency 
aberration that significantly complicates the system of beam 
position control. 



Unlike the FS-BF space processing sub-system, which 
implements the standard narrowband routine of complex- 
valued matrix-matrix multiplication, the new full-band SF- 



BF space processor is designed to operate in two following 
sub-stages: 

1. Wideband Beamsteering realized by the A^-element 
banks of frequency independent digital delays which align 
the input signal samples to focus the w-th SF-BF arm on the 
corresponding direction [7], 

2. Staggered Narrowband Beamforming implemented by 
the Synchronous Commutators (SC) and corresponding 
adders to perform time-distributed narrowband spatial 
integration. 

Let us take a more detailed look at the both techniques of 
the new all-digital wideband beamforming. 

Wideband Digital Time-Aligning 

There three techniques of precise digital time-aligning could 
be considered eligible for use in the SKA aperture array 
processor: 
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Figure 6 - Linear interpolation based fractional delaying 



1. Silicon gate based keyed delay banks, 

2. FIR-filtering based, HW -level realizations consisted of 
FIFO structures for integer (coarse-grained) delaying and 
FIR filters for fractional (fine-grained) delaying [8], 

3. Linear interpolation, SW-level implementations that 
imply re-indexing of data samples for integer delaying and 
time -domain linear interpolation of successive samples for 
fractional delaying. 

The latter method looks like the most attractive to be 
implemented on massively multi-core parallel fixed-point 
(integer) processors. 

The leading idea of the linear interpolation based digital 
delaying is very simple. Generally, a given time delay, ris 

t= Ti+ t f , T/^nTs, n = 1,2, ... (5) 

where Ts is the sample time discrete, Tj and r F is the integer 
and fractional part, respectively. As mentioned above, the 
integer delay can easily be implemented by means of re- 
indexation or using a conventional shift register. 

Let us consider a scenario when a reference wavefront 
crosses the time axis of a given space channel between two 
successive signal samples as shown in Figure 6. 

From the picture is evident that in place of fractional 
delaying we can use the sample that linearly interpolated 
between the early and late samples using the following 
simple expression 

S T = aS L + (1 - a)S E (6) 

where 



Following this "substitution" technique, we can reconstruct 
the true wavefront position for any given boresight. 

Staggered Narrowband Beamforming 

The second stage of spatial beamforming in a single-beam 
SF-beamforming arm, or beam-shaping is implemented by 
Na SC's and N F adders. 

The adders in the beamformer are designed to be without 
precision losses, i.e. the data dynamic range of partial two- 
item sums gradually extends from the front -end to the back- 
end until it becomes sufficient, e.g. 8- to 12-bit. 

The digital FIR-filters and corresponding FFT units aimed 
to finish polyphase frequency demultiplexing in the SF- 
beamformer are identical to those in the FS-beamforming 
system. 

Let us count the total number of real-valued multiplications 
the SF-BF consumes to form N B spatial beams for each of 
N F frequency slots. 

If the SF-beamformer uses the linear interpolation method 
in time aligning the input samples then the number of real- 
valued multiplications at the first "space" stage of 
processing is about 

M SFI = 0{N A N B } (8) 

At the second "frequency" stage of processing, the FS- 
beamformer should perform about 

M SF2 = 0{2N B N F \og 2 N F } (9) 

real multiplications. The total number of operations is 



a = (r s - r 7 )/zj= t 2 /t s 



(7) 



M SF = M SF1 + M SF2 = O {N B (N A + 2N F log 2 JV» } (10) 



4. Conclusions 
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We consider the SKA station beamformer as a wideband 
spatial filter to be matched to any kind of expected radio- 
astronomy signals through maximizing the sample SNR at 
the filter outputs at each time discrete. 

The structure of the new beamformer is built following the 
rule of 'spatial integration first' to avoid processing losses 
in the sample SNR's for all kinds of cosmic signals in the 
environment of short bit-grid size digital representation of 
input data streams. 

Beamsteering is implemented by making use of frequency 
independent high-precision digital time-aligning of antenna 
signals that reduces the complexity of the beamformer and 
keeps the unified time scale unbroken. 

Numerous simulations have shown that as expected the SF- 
beamformer provides the same detection performance of 
weak continuous harmonic signals in additive white 
Gaussian noise (AWGN) as the BF-beamformer does. 

Meanwhile, we expect the SF-beamformer configuration to 
outperform the FS-BF structure in terms of detection 
sensitivity for wideband and full-band signals of interest. 

The newly designed beamforming configuration offers 
much better advantages from the standpoint of effective 
implementation on dedicated ASIC's, or on massively 
multi-core parallel integer ASIP's. 

Using (1) in (4) and (10), we have the following result 



TABLE I 





FS-beamformer 


SF-beamformer 


Real-valued 
multiplications 


-3.6-10 11 


-1.210 s 



Providing the same signal detection capability, the new 
beamforming configuration outperforms by far the known 
structure in terms of computational complexity and hence in 
terms of cost, reliability and power consumption. 
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